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Abstract 

A new site percolation model, directed spiral percolation (DSP), under both directional and 
rotational (spiral) constraints is studied numerically on the square lattice. The critical percolation 
threshold p c as 0.655 is found between the directed and spiral percolation thresholds. Infinite perco- 
lation clusters are fractals of dimension df ~ 1.733. The clusters generated are anisotropic. Due to 
the rotational constraint, the cluster growth is deviated from that expected due to the directional 
constraint. Connectivity lengths, one along the elongation of the cluster and the other perpendicu- 
lar to it, diverge as p — > p c with different critical exponents. The clusters are less anisotropic than 
the directed percolation clusters. Different moments of the cluster size distribution P s (p) show 
power law behavior with \p — p c \ in the critical regime with appropriate critical exponents. The 
values of the critical exponents are estimated and found to be very different from those obtained 
in other percolation models. The proposed DSP model thus belongs to a new universality class. 
A scaling theory has been developed for the cluster related quantities. The critical exponents 
satisfy the scaling relations including the hyperscaling which is violated in directed percolation. A 
reasonable data collapse is observed in favour of the assumed scaling function form of P s (p) . The 
results obtained are in good agreement with other model calculations. 



I. INTRODUCTION 



Recently, there is strong interest in studying properties of electro-rheological and 
magneto-rheological fluids pfl, magnetic semiconductors 2], and composite materials jsj 
because of their industrial applications. A new site percolation model, directed spiral perco- 
lation (DSP), proposed here could be applied in general to study their rheological, electrical 
and magnetic properties. The DSP model is constructed imposing both directional and ro- 
tational constraints on the ordinary percolation (OP) model |4j. The directional constraint 
is in a fixed direction in space and the empty sites in that direction are accessible to occu- 
pation. Due to the rotational constraint the sites in the forward direction or in a rotational 
direction, say clockwise, are accessible to occupation. The direction of the rotational con- 
straint is not fixed in space and it depends on the direction from which the present site is 
occupied. For charged particles, the directional constraint in the model could arise from 
an electric field along a particular direction in the plane of the lattice and the rotational 
constraint may arise from a magnetic field applied perpendicular to the plane. The model 
will also be applicable to the physical situations corresponding to the presence of other kinds 
of directional and rotational constraints. 

The effect of two different external constraints, directional constraint and rotational (spi- 
ral) constraint, on the ordinary percolation model have been studied independently. The 
corresponding models are known as directed percolation (DP) model ^] and spiral percola- 
tion (SP) model ,6]. The DP model has applications in self-organized criticality [7], reaction 
diffusion systems a], nonlinear random resistor networks etc. The SP model has been 
applied in studying spiral forest fire lioj . pinning of interfaces and diffusion under ro- 



tational bias in disordered systems |l2j. It is observed that in the presence of an external 
constraint the critical properties of the system as well as the universality class of the OP 
model are changed The critical exponents associated with different cluster related 

quantities in the DP and SP models are not only different but also different from those in 
the case of the OP model. The DP and SP models then belong to two different universality 
classes than that of the OP model. 

In this paper, the proposed DSP model is studied numerically on the square lattice in two 
dimensions {2D). The external constraints in the model determine the nearest empty sites 
available for occupation. A single cluster growth Monte Carlo (MC) algorithm is developed 



for this model under both the external constraints. Critical percolation probability p c , at 
which a spanning (or infinite) cluster appears for the first time, is determined. The spanning 
clusters are found fractals. The clusters generated in the OP and the SP models are isotropic. 
In DP, the clusters are anisotropic and elongated along the directional constraint. The 
clusters generated in the case of the DSP model are anisotropic but they grow in a direction 
different from the directional constraint due to the presence of the rotational constraint. A 
new type of cluster is thus generated in this model. The connectivity length exponents v\\ 
and u± are estimated. u\\ is approximately equal to the connectivity exponent v of OP but 
both v\\ and v± are different from those obtained in DP. The clusters are less anisotropic 
than DP clusters. Different moments of the cluster size distribution P s (p) become singular as 
p —>■ p c with their respective critical exponents. The values of the critical exponents obtained 
here are different from those obtained in other percolation models like OP, DP, and SP. The 
DSP model then belongs to a new universality class. A scaling theory is developed. The 
critical exponents satisfy the scaling relations within error bars including the hyperscaling 
which is violated in DP. A reasonable data collapse is observed in support of the assumed 
scaling function form of P s {p)- The results obtained in this model are in good agreement 
with the results obtained in the study of magnetoresistance in a model of 3-constituents 
composite material 



II. THE DSP MODEL 



A square lattice of size L x L is considered. A directional constraint from left to right and a 
clockwise rotational constraint are imposed on the system. Due to the directional constraint 
any empty site on the right of an occupied site could be occupied in the next Monte Carlo 
(MC) time step. Due to the rotational constraint the empty sites in the forward direction or 
in the clockwise direction can be occupied. To generate clusters under these two constraints 
a single cluster growth algorithm is developed following the original algorithm of Leath^^. 
In this algorithm, the central site of the lattice is occupied with unit probability. All four 
nearest neighbors to the central site can be occupied with equal probability p in the first 
time step. As soon as a site is occupied, the direction from which it is occupied is assigned 
to it. In the next MC time step, two empty sites due to the rotational constraint and the site 
on the right due to the directional constraint will be eligible for occupation for each occupied 



site in the previous time step. This is illustrated in Fig. The directional constraint is 
represented by two long arrows from left to right. The presence of the rotational constraint 
is shown by the encircled dots. The black circles represent the occupied sites and the open 
circles represent the empty sites. The empty nearest neighbors of the central occupied site 
will be selected for occupation in this step. The direction from which the central site is 
occupied is represented by a short thick arrow. The dotted arrow indicates the eligible 
empty site for occupation due to directional constraint and the thin arrows indicate the 
eligible empty sites for occupation due to the rotational constraint. Since the directional 
constraint is to the right, site 3, an empty site on the right of the occupied site, is always 
eligible for occupation. The sites accessible to occupation due to the rotational constraint 
will be identified now. The rotational constraint acts in the forward or in the clockwise 
direction and depends on the direction of approach to the present occupied site. In Fig. [T] 
(a), the central site is occupied from site 1 on the left. Thus, site 3 in the forward direction 
and site 4 in the clockwise direction are the eligible sites for occupation due to the rotational 
constraint. In this situation, sites 3 and 4 are then the only eligible sites for occupation in 
the next time step due to both the constraints. In (6), the central site is occupied from site 
2 on the top and thus, sites 4 and 1 are the eligible sites for occupation due to the rotational 
constraint. The available sites for occupation are 1, 3 and 4 due to both the constraints in 
this case. In (c), the eligible sites for occupation are 1 and 2 due to the rotational constraint 
and site 3 due to the directional constraint. Note that, site 3 is already an occupied site 
and could be reoccupied from a different direction. A site is forbidden for occupation from 
the same direction. On the square lattice, a site then could be occupied at most 4 times 
from 4 possible directions. This is unlike the case of the ordinary and directed percolations 
where a site is occupied only once. In (d), sites 2 and 3 are eligible for occupation due to 
both constraints. It can be seen that the direction of the rotational constraint is not fixed 
in space and depends on the previous time step whereas that of the directional constraint 
remains fixed in space. In that sense, the directional constraint is a global constraint and 
the rotational constraint is a local constraint in the model. After selecting the eligible sites 
for occupation, they are occupied with probability p. The coordinate of an occupied site in 
a cluster is denoted by (x,y). Periodic boundary conditions are applied in both directions 
and the coordinates of the occupied sites are adjusted accordingly whenever the boundary is 
crossed. At each time step the span of the cluster in the x and y directions L x = x max — x m i n 



and L y = y max — Umin are determined. If L x or L y > L, the system size, then the cluster 
is considered to be a spanning cluster. The critical percolation probability p c is defined as 
below which there is no spanning cluster and at p = p c a spanning cluster appears for the 
first time in the system. 

A typical spanning or infinite cluster generated on a 128 x 128 square lattice at p = 0.655 
is shown in FigEl The black dots represent the occupied sites. The cross near to the 
upper left corner is the origin of the cluster (it was the central site of the lattice). The 
thick arrows from left to right at the top and bottom represent the directional constraint. 
The encircled points in the upper right and lower left corners represent the existence of the 
clockwise rotational constraint. Notice that all the dangling ends are clockwisely turned 
as it is expected. The cluster is highly rarefied. Holes of all possible sizes are there. The 
elongation of the cluster is almost along the left upper to the right lower diagonal of the 
lattice and not along the directional constraint applied in the x direction. The geometry 
of the infinite clusters in other percolation models are the following. In OP, the clusters 
are isotropic. .„ DP, the anisotropic cluster, are elongated along the applied neldQ. In 
SP, the clusters are compact and isotropic 17]. Not only the geometry but also the internal 
structure, highly rarefied with curly dangling ends, of the cluster is different from that of 
the infinite clusters in other percolation models. A new type of anisotropic cluster is thus 
generated in the DSP model. 



III. SCALING THEORY 



The cluster related quantities and their singularities at p = p c are defined here. Scaling 
relations among the critical exponents describing the singularities of the cluster related 
quantities will also be established. The cluster size distribution is defined as 

Psip) = £~ (1) 

where N s is the number of s-sited clusters in a total of N to t clusters generated. In the single 
cluster growth method, the origin is occupied with unit probability. The scaling function 
form of the cluster size distribution is then assumed to be 

Psip) = s- T+1 f[s^p - p c )] (2) 

where r and a are two exponents. 



The order parameter of the percolation transition is the probability that a site belongs 
to a spanning (or infinite) cluster and here it is defined as P^ = p — pJ2' s Ps{p)- The sum is 
over all the finite clusters and it is indicated by a prime. The leading singularity of P^ will 
be governed by J2' s Ps{p), the zeroth moment of cluster size distribution. P^ goes to zero as 
p — > p c from the above with a critical exponent (3. The exponent (3 is defined as 

Poo~(p-PcY- (3) 

An important cluster related quantity is the average cluster size \. In the single cluster 
growth method it is the first moment of the cluster size distribution P s (p) and defined as 
X = J2's s Ps{p) where the sum is over the finite clusters. Two next higher moments xi an d 
Xi are also defined as xi — Yl s s2 Ps{p) and xi — J2' s s 3 P s (p)- In the present context these 
higher moments have no physical meaning but might be useful in particular situations. The 
moments %■> Xu an d X2 diverge with their respective critical exponent 7, 5, and 77 at p = p c . 
The critical exponents 7, 5, and rj are defined as 

x~b-pd~ 7 , xi~\p-Pc\~ 5 , & X2 ~ \p - Pc\' V ■ (4) 

Since the cluster related quantities are just different moments of the cluster size distribution 
function P s (p) then the critical exponents associated with them are not all independent. 
All the critical exponents could be expressed in terms of the exponents r and a needed to 
describe P s (p) (EqJ^j). It could be shown that the kth moment of P s (p) become singular as 

Ks k Ps(p)~(p- Pc r (k - T+2)/a - (5) 

Putting the value of k, the order of the moment, in Eq|S] one could obtain the following 
scaling relations 

(3 = (r - 2)1 a, 7 = (3 - r)/a, 5 = (4 - r)/a, & V = (5 - r)/a. (6) 

Eliminating r and a, scaling relations between f3, 7, 5, and rj could be obtained as 

5 = (3 + 2-f & 77 = 25-7. (7) 

Since the clusters generated here are anisotropic in nature, two lengths, £y and are 
needed to describe the connectivity of the occupied sites. £11 is the connectivity length 
along the elongation of the cluster and £j_ is the connectivity length along the perpendicular 



direction to the elongation. To measure £y and £j_ the moment of inertia tensor T, a 2 x 2 
matrix here, is calculated. For a s-sited cluster, the xy component of the tensor is given by 

s 

T xy = ^Zi.x i -x () ){y i -y () ) (8) 
e=i 

where xt and yi are the x and y coordinates of the £th site and (x ,y ) is the coordinate 
of the center of mass of the cluster. The radii of gyration R\\(s) and R±(s) with respect to 
two principal axes could be obtained as Rf\(s) = Ai/s and R±(s) = A2/S where Ai is the 
largest eigenvalue and A 2 is the smallest eigenvalue of the 2x2 moment of inertia matrix 
T. R± is about the axis passing through (xo,yo) and along the elongation of the cluster 
and R\\ is about the axis perpendicular to the elongation and passing through (x ,y ). The 
connectivity lengths now can be determined as 

The correlation lengths £y and £j_ diverge with two different critical exponents v\\ and v>± as 
p — > p c . The critical exponents u\\ and u± are defined as 

Z\\~\P-Pc\- U " & S±~\P-Pc\- V± . (10) 

The cluster mass is given by the number of sites s in the cluster and is expected to scale as 
s R\\R± f ^ at p — Pc, and it should go as s ^ _R||i?| | d " 1 ' ) above p c , where d is the spatial 
dimension of the lattice and df is the fractal dimension of the infinite clusters generated on 
the same lattice. The percolation probability is the ratio of the number of sites on the 
infinite cluster to the total number of sites, 



R\\Rf f - l) 

tX\\IX 1 



(11) 



for < £|| and < Assuming R\\ ~ ^ and Rj_ ~ two hyperscaling relations could 
be obtained as 

VL{d , _ 1) + „ || = I fc (d-lK-^^-l). (12) 

Eliminating d/ from the above relations another scaling relation could be obtained as 

(d - l)u± + v\\ = 7 + 2(3 = 25 - 37 = 3r] - A5. (13) 

In the following, the values of the critical exponents will be determined and the scaling 
relations will be verified. 



IV. RESULTS AND DISCUSSIONS 



Simulation is performed on a square lattice of size 2 10 x 2 10 . The critical probability p c 
at which a spanning cluster appears for the first time in the system is determined first. The 
probability to have a spanning cluster is given by 

Ps P = ^ L = 1-H'Ps{p) (14) 

^tot s 

where n sp is the number of spanning clusters out of total N tot = 10 4 number of clusters 
generated. In FigJ31 P sp is plotted against the probability of occupation p. Note that, P sp is 
not going to zero sharply at a particular value of p. This is due to the finite size of the lattice 
chosen here. The critical probability p c is then determined from the maximum slope of the 
curve P sp versus p. In the inset of FigElthe slope dP sp / dp is plotted against p. It can be seen 
that the threshold p c is at 0.655±0.001 corresponding to the maximum slope. The derivative 
is calculated using the central difference method for the data points collected in an interval 
of 0.001. Note that, the value of p c obtained here for the DSP model is slightly above the 
value of p c of the directed percolation p c (DP) « 0.6445 [3] and below the spiral percolation 
threshold p c (SP) ~ 0.712 17]. It is expected. Because, the rotational constraint tries to 
make the clusters compact by reoccupying the occupied sites of the directed clusters and, 
at the same time, the directional constraint tries to elongate the compact spiral clusters. 

The infinite cluster, shown in Fig. El generated on a 128 x 128 square lattice has holes 
of almost all possible sizes. It seems that the infinite clusters are self-similar and fractals. 
The fractal dimension df of the infinite clusters generated on the original 2 10 x 2 10 lattice is 
determined by the box counting method. The number of boxes iV^(e) is expected to grow 
with the box size e as Nb(c) ~ t df where df is the fractal dimension. In Fig. Ns(e) 
is plotted against the box size e. The data are averaged over 512 samples. A reasonably 
good straight line is obtained in the log — log scale. The fractal dimension is found df = 
1.733±0.005. The error is due to the least square fitting of the data points taking into account 
the statistical error of each point. The fractal dimension df ~ 1.733 (~ 12/7) obtained here 
is the smallest among the fractal dimensions obtained in other percolation models. The 
values of df in other percolation models are: df{OP) = 91/48 18], df(DP) m 1.765 [19], 



otner per c 
1.957 |lj. 



and df(SP) ~ 1.957 |17(. Also notice that the fractal dimension df(DSP) obtained here 
is little higher than the fractal dimension 1.64 of ordinary lattice animalsj^cj, large OP 



clusters below p c . The infinite cluster generated in the DSP model is then the most rarefied 
one among the infinite clusters obtained in all four models. 

Next, the values of the critical exponents 7, 5, and 77 are estimated. The average cluster 
size x an d two other higher moments Xi an d X2 are measured generating 10 4 finite clusters 
below p c for six different p values. In Fig. |3J Xi Xii an d X2 are plotted against \p — p c \. The 
circles represent Xi the squares represent Xi an d the triangles represent X2- The values of 
the exponents obtained are 7 = 1.85 ± 0.01, 5 = 4.01 ± 0.04, and 77 = 6.21 ± 0.08. The 
errors quoted here are the standard least square fit error taking into account the statistical 
error of each single data point. The values of the exponents 7, 5, and rj are also determined 
by the same Monte Carlo technique on the square lattices of three different smaller sizes 
2 9 x 2 9 , 2 8 x 2 s , and 2 7 x 2 7 to check the finite system size effects on the data. The values 
of the exponents 7 (circles), 5 (squares), and rj (triangles) for different system sizes are 
plotted against the system size 1/L in the inset of Fig. El It is then extrapolated upto 
L — > 00, the infinite system size. The extrapolated values of the exponents are marked 
by crosses and the numerical values obtained are 7 1.85, 5 m 4.01, and rj « 6.21, the 
same as that of the system size 2 10 x 2 10 . These values of exponents are now used to 
verify the scaling relations obtained in section III. One of the scaling relations 77 = 25 — 7 
obtained in EqJT] is checked now. The value of 25 — 7 = 6.17 is very close to the value 
of the exponent rj = 6.21. The scaling relation 77 = 25 — 7 is then satisfied within error 
bars. Using the scaling relation 5 = (3 + 27 (EqEJ), the value of the exponent (3 is obtained 
as (3 = 0.31 ± 0.06. The error has propagated from the errors of 7 and 5. Independent 
estimation of the exponent (3 becomes difficult because of the presence of curvature in the 
plot of versus (p — p c ) in the log-log scale for p > p c . This may be due to dominating 
corrections to scaling to the leading singularity of P^. The values of the exponents r and a 
can also be estimated using the scaling relations in Eq. El Three different values of r could 
be obtained as n = (35 - 4r/)/(S - 7) = 2.14 ± 0.18, r 2 = (4nj - 55)/(r] - 5) = 2.18 ± 0.28, 
and T3 = (3?7 — 5 r y)/(j] — 7) = 2.15 ± 0.11. The estimate of r can be taken as the average of 
Ti, t 2 and r 3 and it is given by r = 2.16 ± 0.20. Similarly a = 0.459 ± 0.015 is determined 
from a x = 1/(5 - 7) = 0.463 ± 0.011, a 2 = 1/(77 - 5) = 0.455 ± 025, and a 3 = 2/(77 - 7) = 
0.459 ± 0.009. The errors quoted here are the propagation errors. A comparison between 
the values of the exponents obtained in different percolation models is made in Table It 
can be seen that the values of the exponents obtained in this model are very different from 



those obtained in other percolation models like OP, DP, and SP. The magnitude of (3 is the 
largest and the magnitude of 7 is the smallest among the four models. The DSP model then 
belongs to a new universality class. The values of the exponents obtained here could be 
approximated to the nearest rational fractions as (3 ~ 1/3, 7 ~ 11/6, 5 ~ 24/6, rj ~ 37/6, 
t Pd 28/13 and a ~ 6/13. Surprisingly they satisfy all the scaling relations in Eq. Inland Eq. 
[7| exactly. 

The connectivity lengths, £y and for the system size 2 10 x 2 10 are plotted against 
\p — p c \ in Fig® The squares represent £11 and the circles represent The corresponding 
exponents u» and u± are obtained as u» = 1.33 ± 0.01 and u± = 1.12 ± 0.03. The errors 
quoted here are the least square fit errors. In the inset of Fig. v\\ (squares) and u± 
(circles) are plotted against different system sizes 1/L. The exponent values are extrapolated 
upto L — > 00 and the extrapolated values are v\\ m 1.33 and u± « 1.12, the same as 
obtained for the system size 2 10 x 2 10 . There are few things to notice. First, the value 
of is almost equal to the connectivity exponent of ordinary percolation, v(OP) =4/3 
jisl ]. Recently, a study of magnetoresistance of a three-component composites consisting of 
cylindrical insulator and perfect conductors in a metallic host film is made by Barabash et 
al Q|. Since the Hall effect will generate an electric field with a component perpendicular 
to the plane of the film with an in-plane applied current, their system appears inherently 
three dimensional (3D). However, the electric field perpendicular to the film plane vanishes 
because of the presence of columnar perfect conductor. Thus, their 3D problem reduces 
to that of calculating the effective conductivity of a 2D composites of perfect insulator 
and perfect conductor jjj Q]. The results obtained in the model of composite material by 
Barabash et al [ljj then could be compared with the results of the present DSP model in 
2 dimensions. It is found by Barabash et al Q| that the correlation length exponent is 
4/3 independent of anisotropy. The correlation length exponent quoted there is equivalent 
to the connectivity exponent u» of the DSP model considered here. The results of the 
DSP model is thus in good agreement with the results obtained in the model calculation of 
magnetoresistance of composite materials. Second, the values of u» and u± are different from 
those obtained in the DP model, V\\{DP) = 1.733±0.001 and u L (DP) = 1.0972±0.0006 |2l| 
(see Table H] also). Third, the ratio of the connectivity lengths goes as ~ \p ~ Pc\~ Au 

where Au = u\\ — v±. In DSP, Av is approximately 0.21 whereas in DP, it is approximately 
0.64. This means that clusters in the DSP model are less anisotropic than the clusters in the 



DP model at any p. This is due to the presence of the rotational constraint which tries to 
make the clusters isotropic. Barabash et al ^| assumed that Av = 0. The effective aspect 
ratio of their problem might vary as \p — p c \~ Au with Av = 0.21 as obtained here in the DSP 
model. Fourth, (d — l)v± + v\\ is ~ 2.45 for d = 2 and 25 — 37 ~ 2.47. The hyperscaling 
relation 25 — 37 = {d — + v\\ in Eq. EI] is then satisfied within error bars. Two other 
hyperscaling relations v±(df — 1) + Vn = I /a and (d — l)v± — (3 = v±(df — 1) given in 
El are also satisfied within error bars. In directed percolation hyperscaling is violated 
This is the first anisotropic percolation model where hyperscaling is satisfied. The values of 
the connectivity exponents v» and v± could also be approximated to their nearest rational 
fractions as V\\ ~ 4/3 and v± ph 7/6. They satisfy the hyperscaling relations in Eq. fT2l and 
[EH exactly taking df = 12/7. 

In the above study, it is found that the fractal dimension df and the values of the critical 
exponents (r, a, /3, 7, v\\, vj_, etc.) of DSP model are different from other percolation 
models. As a consequence, the DSP model belongs to a new universality class. This is due 
to the fact that a completely new type of percolation cluster is generated in the DSP model. 
There are three important features of the DSP clusters. The clusters are rarefied, have 
spiraling danging ends, and are anisotropic. As the clusters grow, more and more vacancies 
are generated into the cluster. This is because, the clockwise rotational constraint tries to 
occupy sites away from the directional constraint whereas the directional constraint tries to 
occupy sites along itself in the x-direction. The effects of these two counter acting constraints 
are the following: the spiral clusters (clusters with spiraling dangling ends) become rarefied, 
become wider as it is away from the origin (see FigJ3) , and elongated along a clockwisely 
rotated direction from the original globally fixed directional constraint. The clusters are 
neither directed percolation clusters nor spiral percolation clusters. The combined directed 
and spiral constraints produces highly rarefied anisotropic spiral clusters. Individually, the 
effect of these three features, anisotrpoy [16], spiraling [17], and volume fraction 23], on 
the ordinary percolation clusters have been studied. Each of them corresponds to different 
critical behaviour and consequently belongs to new universality class. Here in DSP model, 
a new critical behaviour is obtained at the percolation threshold because of the presence of 
all three features in the same cluster. 

Finally, the scaling function form assumed for the cluster size distribution P s (p) = 
s~ T+1 f [s a (p — p c )] is verified. The scaled cluster size distribution P s (p) / P s (p c ) is now plotted 




against the scaled variable s a (p — p c ) in Fig. [7| It is assumed that f [0] is a constant. The 
value of a is taken as 0.459. Distribution of 10 4 clusters over the bins of width 2* _1 — (2* — 1) 
with % = 1 to L 2 (= 2 20 ) is considered for each p value. In Fig. the cluster size s changes 
from 64 to 8192 and (p — p c ) changes form 0.12 to —0.06. Data for different values of s and 
(p ~ Pc) collapse reasonably onto a single curve. This means that the scaling function form 
assumed for the cluster size distribution is appropriate for this model. Also notice that the 
maximum value of P s (p)/P s (p c ) for DSP model (~ 3) is different from that of OP model 
(«4.5) Q. 



V. CONCLUSION 



The new directed spiral site percolation, DSP model, under both directional and ro- 
tational constraints belongs to a new universality class. The combined directed and spiral 
constraints produces a new type of percolation clusters which are highly rarefied, anisotropic, 
and spiral in nature. The critical properties of the cluster related quantities in this model 
are very different from the other percolation models like OP, DP, and SP. The critical expo- 
nent 7 of the average cluster size % is the smallest and the exponent f3 of is the largest 
among the four models. The fractal dimension df « 12/7 is the smallest among the four 
percolation models and thus the infinite clusters are highly rarefied here. Since the clusters 
generated in this model are anisotropic, two connectivity lengths, £y and are defined to 
describe the scaling behaviour of the clusters connectivity as p — > p c . The exponents u» 
and u±, associated with £m and are estimated. The finite size effect on the exponent 
values is also checked making simulations on different system sizes and the values of the 
exponents are extrapolated to the infinite network. It is found that uu is approximately 
equal to the connectivity exponent v = 4/3 of the OP model which is also in good agree- 



ment with the correlation length exponent v obtained by Barabash et al 14| in the study of 
magnetoresistance of a model composite material. Both uu and u± are different from those 
obtained in the DP model. The order of anisotropy £y/£± is higher in the DP clusters than 
that of the clusters in the DSP model at any p. The critical exponents of the DSP model 
satisfy the scaling relations within error bars including the hyperscaling which is violated 
in the DP model. This is the first anisotropic percolation model where the hyperscaling is 
satisfied. The assumed scaling function form P s (p) = s~ T+1 f [s a (p — p c )] of the cluster size 



distribution P s (p) is verified through data collapse. The values of all the critical exponents 
are suggested in terms of rational fractions and it is observed that those rational fractions 
satisfy the scaling relations exactly. The proposed rational fractions for the values of the 
critical exponents might be verified through an exact solution of the model. A difficulty 
in solving the model exactly is in the reoccupation of the same sites from different direc- 
tions. The model will be applicable to the physical situations where both the directional 
and rotational constraints are present. For example, measurement of magnetoresistance in 
composite materials, magnetic semiconductors, and super-ionic conductors could be studied 
using this model if it is extended to 3 dimensions. The model will also be applicable in 
studying the rheological properties of electro-rheological and magneto-rheological fluids in 
presence of crossed electric and magnetic fields. 
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Percolation 
Models 
OP 
DP 



7 



a 



18 
19 



21] 



SP 



Q 



DSP 



5/36 


43/18 


187/91 


36/91 


4/3 


91/48 


0.277 


2.2772 


2.108 


0.3915 


v± = 1.0972 


1.765 


±0.002 


±0.0003 


±0.001 


±0.0004 


±0.0006 












= 1.733 












±0.001 




0.048 


2.19 


2.022 


0.447 


1.116 


1.957 


±0.011 


±0.07 


±0.004 


±0.014 


±0.003 


±0.009 


0.31 


1.85 


2.16 


0.459 


v± = 1.12 


1.733 


±0.01 


±0.01 


±0.20 


±0.004 


±0.03 


±0.005 


(1/3) 


(11/6) 


(28/13) 


(6/13) 


(7/6) 


(12/7) 










v\\ = 1.33 












±0.01 












(4/3) 





TABLE I: Comparison of the values of the critical exponents /3, 7, r, a, v, and df in the case 
of ordinary (OP), directed (DP), spiral (SP), and directed spiral (DSP) percolation on the square 
lattice. The values of the critical exponents of the DSP model are different from the other models. 
The values within parenthesis are the nearest rational fractions of the values of the critical expo- 
nents. These rational fractions satisfy the scaling relations exactly. The DSP model belongs to a 
new universality class. 
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(d) 



FIG. 1: Selection of empty nearest neighbors of an occupied site for occupation. Black circles 
are the occupied sites and open circles are the empty sites. Two thick long arrows from left to 
right represent the directional constraint. The presence of clockwise rotational constraint is shown 
by the encircled dots. The eligible empty nearest neighbors of the central occupied site will be 
selected here for occupation. The direction from which the central site is occupied is indicated by 
a short thick arrow. The dotted arrow indicates the sites allowed by the directional constraint and 
the thin arrows indicate the sites allowed by the rotational constraint. In (a), the central site is 
occupied from the left and the sites 3 and 4 are accessible to occupation. In (b), the central site is 
approached from the top and sites 1, 3 and 4 could be occupied. In (c), the central site is occupied 
from the right and 1, 2, and 3 are probable for occupation. Notice that the site 3, already an 
occupied site, could be reoccupied from a different direction. In (d), the central site is approached 
from the below and sites 2 and 3 are the probable sites for occupation. 
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FIG. 2: An infinite cluster on a 128 x 128 square lattice at p = 0.655 is shown. The black dots are 
the occupied sites. The cross on the upper left corner is the origin of the cluster. The thick arrows 
on the top and bottom from left to right represent directional constraint. The presence of the 
clockwise rotational constraint is shown by the encircled dots. The cluster is highly rarefied and 
has holes of almost all possible sizes. The elongation of the cluster is along the upper left to the 
lower right diagonal and not along the directional constraint. The dangling ends are clockwisely 
rotated. 
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FIG. 3: Plot of spanning probability P sp versus p (circles). In the inset the slope dP sp /dp is plotted 
against p (squares). The critical probability p c is determined from the maximum slope and it is 
found p c = 0.655 ± 0.001 as indicated by an arrow. 
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FIG. 5: Plot of the first, second and third moments Xi Xii an d X2 of the cluster size distribution 
versus \p — p c \. Different symbols are: circles for x, squares for xi, and triangles for X2- The 
corresponding critical exponents are found as 7 = 1.85 ±0.01, 5 = 4.01 ±0.04, and rf = 6.21 ±0.08. 
In the inset, the values of the exponents 7(0) > <5( a ) an d 77(A) are plotted against the system size 
1/L. Extrapolating to the infinite system size (1/L = 0), the values of the exponents obtained are 
7 1.85, 5 w 4.01, and r) 6.21. 
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FIG. 6: The connectivity lengths, £y and are plotted against \p — p c \. The circles represent 
and the squares represent £y. The critical exponents are found as i/y = 1.33 ± 0.01 and u± = 
1.12 ± 0.03. In the inset, the values of and ^±(0) are plotted against the system size 

1/L. The values of the exponents are extrapolated upto L — > oo. The extrapolated values of the 
exponents are marked by crosses and they are vu ~ 1.33 and ~ 1.12. 




FIG. 7: Plot of the scaled cluster size distribution P s (p) / P s (p c ) versus the scaled variable s a (p— p c ) 
for different values of p with a = 0.459. The cluster size s changes from 64 to 8192. The data 
plotted correspond to p - p c = 0.12(x), 0.11(A), -0.03(0), -0.04(D), -0.05(O), and -0.06(v)- 
A reasonable data collapse is observed. 



